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Low-lying A baryons with spin 1/2 are analyzed in two-flavor lattice QCD. In order 
to extract two low-lying states for each parity, we construct 2x2 cross correlators from 
flavor SU(3) "octet" and "singlet" baryon operators, and diagonalize them. Two-flavor 
CP-PACS gauge configurations are employed, which are generated with the renormalization- 
group improved gauge action and the C(a)-improved quark action. Simulation are performed 
at three different /3's, /3 = 1.80, 1.95 and 2.10, whose corresponding lattice spacings are a = 
0.2150, 0.1555 and 0.1076 fm. For each cutoff, we adopt four different hopping parameters, 
(re va i,K soa ). The corresponding pion masses range from about 500 MeV to 1.1 GeV. The 
results are extrapolated to the physical quark- mass point. Our results indicate that there 
are two negative-parity A states nearly degenerate at around 1.6 GeV, whereas no state as 
low as /1(1405) is observed. By extracting the flavor components of each state, we find that 
the lowest (lst-excited) negative-parity state is dominated by flavor-singlet (flavor-octet) 
component. 

§1. Introduction 

;1(1405) has been attracting much interest from several view points. A(1405) is 
the lightest negative-parity baryon in spite of the valence strange quark in it. Among 
the J p = l/2~ baryons, 4(1405) is much lighter than the non-strange counterpart 
iV(1535), and isolated from the others. It has no spin-orbit partner in the vicinity, 
taking into account that the lowest spin 3/2~ state is /1(1520). Furthermore, the 
structure of 4(1405) remains mysterious. On one hand, 4(1405) is interpreted as a 
navor-SU(3)-singlet three-quark state in conventional quark models. On the other 
hand, 4(1405) could be interpreted as an antikaon-nucleon KN molecular bound 
state (B.E. ~ 30 MeV). The binding energy of KN implies a strong attraction 
between K and iV which may cause a new type of dense hadronic matter, 

kaonic nuclei or kaonic nuclear matter pOSS We also expect that such KN bound 
states with large binding energies can be regarded as compact 5-quark states. The 
5-quark picture of 4(1405) has advantages that all five quarks can be placed in the 
lowest-lying L = state forming a negative-parity baryon, and also that it requires 
no spin-orbit partner for 4(1405). 

Therefore, the property of yl(1405) can be an important clue to new paradigm 
in hadron physics. We here study properties of yl(1405) using the lattice QCD 
formulation. Though several lattice QCD studies on A(1405) have been performed so 
far mm mm most f them are based on quenched QCD, and few lattice QCD studies 
succeeded in reproducing the mass of 4(1405). Moreover, in many cases, individual 
analyses with "singlet" and "octet" operators were performed, where possible flavor 
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mixings were not taken into account. Little has been discussed on the lattice about 
the possible mixing of flavor-SU(3)-octet and -singlet components induced by the 
symmetry breaking. It is then an intriguing issue to clarify the flavor structures in 
excited-state hadron resonances. 

Several possible reasons for the failure of reproducing /1(1405) in lattice QCD 
were suggested through these studies; lack of meson-baryon components due to 
quenching, exotic (non-3-quark type) structure of 71(1405), or insufficiency of the 
lattice volume in the simulations. Resolving such difficulties requires unquenched 
lattice QCD calculation on a larger lattice volume with varieties of interpolating 
operators. Evaluating contaminations by scattering states induced purely by sea 
quarks would be newly needed. In this paper, we aim at clarifying the properties 
of yl(1405) with two-flavor full lattice QCDpB employing the "octet" and "singlet" 
baryon operators to construct correlation matrices, which enables us to extract the 
low-lying spectrum as well as the mixing between octet and singlet components in 
71(1405). 

§2. Lattice QCD setups 

2.1. Simulation conditions 

We adopt the renormalization-group improved gauge action and the C(a)-improved 
quark action. We adopt three different /3's, /3 = 1.80, 1.95 and 2.10, and corre- 
sponding lattice spacings are a = 0.2150, 0.1555 and 0.1076 fmP We employ four 
different hopping parameters (K V ai>^sea) for each cutoff. Corresponding pion masses 
range approximately from 500MeV to 1.1 GeV at each f3. The details are found in 
Ref.IUD 

2.2. Baryonic operators 

In order to extract the low- lying states in S = —1 and isosinglet channel, we 
construct 2x2 cross correlators from the "singlet" and "octet" operators, 

Vl (x) ee -Le abc 

x [u a {x)[d Tb {x)C lb s c {x)} + d a {x)[s Tb {x)C lb u c {x)} + s a (x)[u Tb {x)C lb d c (x)}) 

x (u a {x)[d Th {x)C lb s c (x)} + d a {x)[s Tb {x)C lb u c {x)} - 2s a (x)[u Tb (x)Cj 5 d c {x)]j 

It is easy to check that ?yi(x) (r]s(x)) belongs to the singlet (octet) irreducible rep- 
resentation of the flavor SU(3). We adopt point-type operators for the sink, and 
extended operators, which are smeared in a gauge-invariant manner, for the source. 
Smearing parameters are chosen so that root-mean-square radius is approximately 
0.5 fm. 
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§3. Lattice QCD Results 



3.1. Hadronic masses 

Fig. [T] (left and middle panels) show the eigen-energies in the positive- and 
negative-parity channels, plotted as functions of the squared pion mass m^, which 
were obtained on 24 3 x 48 lattice. The lattice spacing is a = 0.1076 fm and the 
lattice cut off is a -1 = 1.83 GeV. Filled circles and open squares denote the energies 
of the ground states and the lst-excited states, respectively. The solid curves rep- 
resent quadratic chiral hts as functions of squared pion mass m^. Solid lines at the 
vertical axes indicate the empirical masses of yl(1115) and yl(1600) in the left panel 
(positive-parity states), and those of /t(1405), 71(1670) and vl(1800) in the middle 
panel (negative-parity states). 

The extrapolated values of the positive-parity ground state agree very well with 
the mass of the ground-state yl(1115). On the other hand, the lst-excited state in this 
channel lies much higher than vl(1600), which is the 1st excited state experimentally 
observed so far. The same tendency was reported in Ref. [8]), and the situation 
is similar to the case of the Roper resonance, the non-strange SU(3) partner of 



In the negative-parity channel (middle panel in Fig{TJ, the ground- and the 
lst-excited states always have similar energies at all the k's. These eigen-energies 
have similar quark-mass dependences, and the mass splittings are almost quark mass 
independent. The chirally extrapolated values both lie around the mass of yl(1670) 
rather than yl(1405). Similarly to the previous studies, the mass of yl(1405) is not 
reproduced in our calculation. While in quenched simulations in Refsj6|),[7|) such 
failure was regarded as an evidence of possible meson-baryon molecule components 
in yl(1405), our present unquenched simulation contains effects of dynamical quarks 
and thus should incorporate meson-baryon molecular states. 

We here mention possible contaminations from scattering states. Since our cal- 
culations contain dynamical up and down quarks, scattering states could come into 
the spectra. (NK and E ir thresholds are drawn in Fig. Q] as dashed lines.) The 
energy of the lst-excited state at the lightest u- and d-quark masses is very close to 
the threshold on the lattice, and could be contaminated by such scattering states. 
There are two meson-baryon channels relevant in the present calculation, irU and 
KN . In the physical charge basis, we have 5 different thresholds (2 for KN, 3 for 
irU). In order to distinguish resonance states from all these scattering states, we 
impose the following boundary condition on the quark fields, 



Under such boundary condition for quark fields, a hadronic state (j)^k+n{x) which 
contains 3k + n valence quarks obeys 



A(1600)P) 



ip(x + L) = e^ip(x). 



(3-1) 



^k+n{x + L) 



= es 



(3-2) 



and can have spatial momenta, pi a t 



3gf (m € Z). 



4 



T. T. Takahashi and M. Oka 




Fig. 1. [Left and middle panels] Masses of the lowest two A states plotted as functions of the 
squared pion mass, which were obtained on 24 3 x 48 lattice. The filled circles (open squares) 
denote the masses of lowest (lst-excited) state. Two solid curves represent quadratic functions 
used in the chiral extrapolation, (left): Two dashed lines indicate the ttE and the KN thresh- 
olds in the presence of the relative momentum p = ^£ . Two solid lines on the vertical axes show 
the experimentally observed masses of yl(1115) and yl(1600). (middle): Two dashed lines indi- 
cate the ttE and the KN thresholds with the relative momentum p — 0. Two solid lines on the 
vertical axes show the experimentally observed masses of /1(1405) and 71(1670). [Right panel] 
Lowest two eigen-energies in the (J p , S) = (l/2~ , —1) channel under the normal (periodic) and 
the twisted boundary conditions. The circles (squares) are for the ground (lst-excited) states. 
The lower and upper solid lines respectively represent the threshold energies of Sir and NK 
states evaluated with normal/twisted boundary condition. These data were obtained with the 
largest hopping parameter (smallest quark mass), and the corresponding pion mass is 550 MeV. 

As a result, only states which consist of 3k valence quarks can be zero-momentum 
states. Since other quark combinations inevitably have non-vanishing spatial mo- 
menta, their energies are raised up. As long as we employ three-quark operators for 
baryon creation/annihilation, hadronic states appearing in scattering states should 
contain one or two valence quark(s), since sea-quark pairs themselves cannot carry 
flavors. 

We plot in the right panel of Fig. [T]the lowest two eigen-energies in the ( J p , S) = 
(l/2~ ,— 1) channel under the periodic and the twisted boundary conditions. The 
open and filled circles (squares) are for the ground (lst-excited) states. The lower 
and upper solid lines respectively represent the threshold energies of Utt and NK 
states under normal/twisted boundary condition. These data were obtained with 
the largest hopping parameter (smallest quark mass), and the corresponding pion 
mass is 550 MeV. The threshold energies are raised up in the case of twisted bound- 
ary condition, whereas the lowest two eigen-energies remain unchanged. Thus we 
conclude that the observed states are insensitive to boundary conditions, and con- 
taminations from meson-baryon scattering states are small, which would be due to 
the same reason as the absence of string breaking in heavy-quark potentials from 
Wilson loops. 

3.2. Flavor structures 

The chiral unitary approach^} has suggested that ^4(1405) is not a single pole 
but a superposition of two independent resonance poles. This possible structure 
of yl(1405) is now attracting much interest, and desired to be clarified in a model 
independent manner. We investigate the flavor structures of the ground and the 
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lst-excited states obtained from the cross correlators of two operators. 

In order to clarify the flavor structures of the low- lying states, we evaluate the 
overlaps of the obtained states with the source and sink interpolating fields. We 
evaluate go and g\ defined as 



g = (Oth | r/8 1 vac) /(Oth 1 7/1 1 vac), g 1 = (lth|7/i|vac)/(lth|7/8|vac) 



(3-3) 



Both </q" and vanish when the SU(3)ir symmetry is exact, showing that the ground 
(1st excited) state is purely flavor singlet (octet) in the limit. For the positive-parity 
states, we similarly define 

g^ = (0th 1 7/i | vac) /(0th 1 7/8 1 vac), g± = (lth|7/8|vac)/(lth|7/i|vac) (3-4) 

In this case, we exchange the denominator and the numerator since the flavor as- 
signment in the SU(3)i? limit of the ground- and the lst-excited states are opposite. 
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Fig. 2. [Left and middle panels] gf are plotted as a function of m^.. Solid lines denote quadratic 
functions used in chiral extrapolation, (left): The ratio gf in the positive-parity channel, 
(middle): The ratio g~ in the negative-parity channel. [Right panel] Effective masses of the 
negative-parity A ground state under the normal (periodic) and the twisted boundary conditions. 

We show g^ and gf' as functions of pion-mass squared for each lattice cutoff in 
Fig. (Left and middle panels). As is expected from the symmetry, such mixing 
coefficients cross zero fit the flavor symmetric limit (/^ = Kd = K s ). The data 
show smooth quark-mass dependence toward the chiral limit in all the lattice-cut-off 
cases and the dependences are almost lattice-cut-off independent. We find that the 
magnitudes of operator mixing get larger for larger flavor-symmetry breaking. It is 
interesting that the lst-excited state in positive-parity channel is dominated by the 
flavor-singlet component showing almost no contribution of octet components. 

On the other hand, the mixing of the singlet and octet components is generally 
weaker in the positive parity A's. It is natural because the splitting of the two states 
are large in the SU(3)_f. Nevertheless, it is interesting to see that the ground state, 
corresponding to yl(1115), contains significant (about 10%) flavor singlet component 
in the chiral limit, which was not expected from the simple quark model picture. 



§4. Conclusion and Discussions 



An unquenched lattice QCD study for the low-lying A baryon states in the S = 
— 1, 1 = and J = 1/2 channel has been performed, focusing especially on the masses 
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and flavor structures of the two lowest negative-parity A's. We have constructed 
2x2 cross correlators from the "octet" and "singlet" baryonic interpolating fields, 
and diagonalized them so that we can properly extract the information of the low- 
lying two states. In our measurements, we have adopted the 2-flavor full-QCD gauge 
configurations generated by CP-PACS collaboration with the renormalization-group 
improved gauge action and the C(a)-improved quark action. We have performed 
simulations at three different /3's, f3 = 1.80, 1.95 and 2.10, whose corresponding 
lattice spacings are a = 0.2150, 0.1555 and 0.1076 fm, employing four different 
hopping parameters (K va i,K sea ) for each cutoff, corresponding to the pion masses 
ranging from 500 MeV to 1.1 GeV. The results have been extrapolated to the physical 
quark-mass region. 

It is shown that the mass of the ground-state /1(1115) is reproduced well at 
all three /3's, while the lst-excited positive-parity A state lies much higher than 
the experimentally observed /1(1600). The same tendency was reported in Ref. [8|), 
and also in the case of the Roper resonance, TV* (1440), which is the non-strange 
SU(3) partner of /1(1600)P2J It should also be noted that no evidence of irU or KN 
scattering states are seen in the calculation. 

Our results for negative-parity A states indicate that there are two l/2~ states 
nearly degenerate at around 1.6 — 1.7GeV, while no state as low as /1(1405) is ob- 
served. We have revealed the flavor structures of these states from the lattice data 
for the first time. It is found that the lowest negative-parity A state is dominated 
by flavor-singlet component. The second state, which is less than 100 MeV above 
the ground state, is predominantly flavor octet. Thus we find that the first two 
negative-parity /l's have different flavor structures. 

However, they come significantly heavier than the experimentally observed lowest- 
mass state, yl(1405). A possible straightforward interpretation is that these obtained 
states do not correspond to the physical yl(1405), but describe the higher A reso- 
nances. In fact, the 2nd and 3rd negative-parity states lie at 1670 MeV and 1800 
MeV, both of which are three-star states in the Particle Data Group classification. 
The present lattice-QCD data are consistent with these excited states. In the non- 
relativistic quark model approach, each of these states is classified as a flavor octet 
P-wave baryon. We, however, have shown that the lower state is dominated by a 
flavor-singlet component. The results here predict one flavor-singlet state and one 
flavor-octet state in the vicinity. The reason for missing /1(1405) state will be its 
poor overlap with three-quark operators. In fact, the inclusion of dynamical quarks 
does not strongly enhance the signals of possible meson-baryon scattering states. 
(See the later discussions.) If 71(1405) is predominantly a meson-baryon molecular 
state, such overlaps would be naturally small. 

The other possibility is, of course, that the lowest and the 2nd-lowest states 
describe physical /1(1405) and 71(1670) but the masses have been overestimated. 
Considering that our simulation contains two-flavor dynamical quarks, the failure of 
obtaining a light A state could be attributed either to (1) uncertainty of the chiral 
extrapolation, (2) strange-quark quenching, (3) insufficient lattice volume or (4) lack 
of chiral symmetry. 

If /1(1405) is not a simple 3-quark state but predominantly a meson-baryon 
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molecule state, the dependence of its mass might be more complicated than 
simple polynomials, as in the case of the m^-dependence of the sigma-meson pole 
in TT-ir scattering.^ The strange-quark quenching (2) seems to make the masses 
of octet baryons in positive parity channel slightly (~ 10%) overestimated in the 
present setupsPJ On the other hand, the deficiencies, (3) and (4), may cause the 
lowest state not properly reproduced, supposing that the main component of /1(1405) 
is a meson-baryon molecular state. In order to check whether this conjecture is 
correct, simulations with light dynamical quarks (m^ <C 500 MeV) and larger volume 
(L 3> 2.5 fm) will be required. 

Upon the above conjecture, we can further consider one interesting scenario that 
these states both correspond to the physical yl(1405). Then, the results may support 
its double pole structure proposed by the chiral unitary approach.^ In our results, 
the lowest two states are almost degenerate at all the /3's (lattice spacing) and ft's 
(quark masses). Namely, the obtained two states are the signature of the double- 
pole resonance, but the mass has not yet been reproduced because of the deficiency 
stated above. 

It is also important to confirm that the meson-baryon scattering states are not 
observed. Despite that dynamical up and down quarks are included and the meson- 
baryon thresholds appear around/below the obtained eigen-energies, we have found 
no clear signal of the meson-baryon scattering states. This fact is supported by the 
observation that the effective mass plots under the normal and twisted boundary 
conditions show no prominent difference. We show in the right panel of Figj2] the 
effective mass plots of the ground state negative-parity A under the normal and 
twisted b.c.'s. They are obtained on the 24 3 x 48 lattice with the largest hopping 
parameter. Thus we conclude that no scattering state appears in the present spec- 
trum. Also in Ref. I15p . in which excited-state nucleon spectrum was systematically 
and extensively investigated with two-flavors of dynamical quarks, no clear signal of 
scattering states was found and the importance of multi-quark operators was raised. 

A similar situation can be found in the computation of the Wilson loops, whose 
expectation values give us the potential between a (heavy fundamental) quark and 
an antiquark. In the presence of dynamical quarks, such an interquark potential 
should saturate and flatten at some interquark distance, where the confining string 
is broken and a quark- antiquark pair is created. However, no one has ever observed 
successfully such a "string breaking" effect by means of the Wilson loopP^ One 
possible reason for this phenomenon is that the Wilson loop itself has poor overlaps 
with such broken strings. 

In the A spectrum, we may similarly conjecture that the scattering states or 
broad resonances, which are sensitive to boundary conditions, have little overlap 
with the 3-quark interpolating field operators. More detailed investigation would be 
needed for clarification of scattering states purely induced by dynamical quarks. 

The flavor structures of the A states have been very well clarified using the 
variational method. The octet and the singlet components are mixed when the flavor- 
SU(3) symmetry is broken. Actually, the ground (lst-excited) state is dominated 
by singlet (octet) component, and the contamination by another representation is at 
most 20% (5% when squared) in our present analysis. From these findings, we expect 
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that the flavor-SU(3) symmetry is not largely broken. A similar conclusion was also 
derived for the study of the meson-baryon coupling constants in lattice QCDP^ 

Because the SU(3) breaking effect seems small, the analyses without SU(3) mix- 
ings adopted so far^'^®® make sense to some extent. The mixings, however, get 
larger towards the chiral limit and variational analyses would be essentially needed, 
when we adopt much lighter quarks. We have also found that the meson-baryon 
contents in each state strongly depend on the meson-baryon couplings and their 
signs. (See Ref. fT0|) .) Precise determination of the meson-baryon contents will re- 
quire reliable determination of the couplings up to signs, for which further lattice 
QCD calculations may be helpful. 
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